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Dipolar particles in an elongated trap are expected to undergo a quantum phase transition from 
a linear to a zigzag structure with decreasing transverse confinement. We derive the low energy 
effective theory of the transition showing that in presence of quantum fluctuations the zigzag phase 
can be characterized by a long ranged string order, while the local Ising correlations decay as a 
power law. This is also confirmed using density matrix renormalization group (DMRG) calculations 
on a microscopic model. The non local order in the bulk gives rise to zero energy states localized at 
the interface between the ordered and disordered phases. Such an interface naturally arises when 
the particles are subject to a weak harmonic confinement along the tube axis. We compute the 
signature of the edge states in the single particle tunneling spectra pointing to differences between 
a system with Bosonic versus fermionic particles. Finally we asses the magnitude of the relevant 
quantum fluctuations in realistic systems of dipolar particles, including ultracold polar molecules as 
well as alkali atoms weakly dressed by a Rydberg excitation. 



I. INTRODUCTION 

The realization of ultracold dipolar gases opens new di- 
rections for investigation of quantum many body physics. 
Relevant systems currently under investigation include 
degenerate gases of atoms with large magnetic dipole mo- 
ments, such as Cr^ or of heteronuclear molecules with a 
large permanent electric dipole.^ Another promising pro- 
posal is to use degenerate Alkali atoms which are weakly 
dressed with a Rydberg excitation by optical pumping.^ 

The long range dipolar interactions in these systems 
can be strong enough to drive interesting structural phase 
transitions. But, not so strong as to make the kinetic en- 
ergy negligible. The balance between kinetic and interac- 
tion terms results in strong quantum fluctuations, which 
provide a fertile ground to formation of novel phases.^'^ 

In this paper we shall specifically consider the zigzag 
instability of a chain of repulsive particles in an elon- 
gated trap, shown schematically in Fig. 1. The dipoles 
are assumed to be polarized by an external electric field 
perpendicular to the trap axis. When the transverse con- 
finement is lowered below a critical value, the repulsive 
interactions between dipoles overcomes the confinement, 
leading to a staggered distortion of the chain. A classi- 
cal analysis along these lines^ explains the zigzag distor- 
tion observed in chains of trapped ions.^ In the classical 
description of either the coulomb or dipolar crystal, the 
zigzag distortion is associated with breaking of the Z2 re- 
flection symmetry about the midplane of the trap. How- 
ever, because the dipolar interactions are much weaker, 
the effect of quantum fluctuations in the particle posi- 
tions is greatly enhanced compared to the ions. We shall 
see that this leads to an interesting change of the zigzag 
phase and phase transition. 

One effect of quantum fluctuations is to allow single 
particle tunneling between an up and down displacement 
in the zigzag. It is convenient to think of distortion of 
a particle, up or down as an Ising spin variable. The 
tunneling then acts to disorder the zigzag in the same 
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FIG. 1: (a) One of the two broken symmetry zigzag states of 
the classical system, (b) Longitudinal quantum fluctuations 
affect a continuous distortion from one zigzag configuration 
to the other, which restores the Z2 symmetry of the ground 
state in the bulk 



way that a transverse field acts to disorder a one dimen- 
sional Ising model inducing transitions between the up 
and down state of the spin. This class of fluctuations and 
a mapping to an Ising model were recently discussed. ^'^ 

Another mode of quantum fluctuation, which has not 
been considered in this context so far, is the axial motion 
of the dipolar particles. In one dimension these longitu- 
dinal fluctuations lead to disordering of the crystal and 
at the same time they restore the Z2 symmetry of the 
staggered configuration. What is then the fate of the 
phase transition and of the broken symmetry phase in 
the presence of the fluctuations? 

To answer this question we derive a long wavelength 
theory, which describes both the gapless longitudinal 
fluctuations and the transverse distortion. Let us briefly 
summarize the main results of this analysis. First, we 
show that the zigzag state remains a distinct phase 
in presence of the fluctuations, however the zigzag or- 
der becomes non-local and is described by a string 
order-parameter. The longitudinal density fluctuations 
(breathing modes of the zigzag) remain gapless in the 
zigzag phase, while single particle excitations induce do- 
main walls in the string order and are therefore gapped. 
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We further show that the presence of string order in the 
bulk imphes the existence of a zero energy single-particle 
state localized at each edge. These states are represented 
by localized Majorana modes in the low energy theory. 
We compute the local tunneling spectra into the edge and 
identify the signatures of the localized zero energy states. 
Finally we point to interesting differences in these signa- 
tures depending on whether the dipolar particles are of 
a Bosonic or a fermionic species. 

Before proceeding we note previous analysis done in 
the context of electronic quantum wires, ^^^^ which also 
identify a quantum phase transition, of the same univer- 
sality class as we consider here, into a zigzag chain. While 
the non-local nature of the order parameter is implicit in 
these studies, this non-locality has not been emphasized, 
and in particular its implications on the edge structure 
have not been discussed. We also note earlier studies on 
Bosonic^^ and Fermionic^^ ladders, where a state simi- 
lar to the zigzag is identified as an out of phase charge 
density wave between the two legs of the ladder. 

The rest of the paper is organized as follows. In sec- 
tion II we identify the important scales in the problem of 
dipolar particles in an elongated trap and derive an effec- 
tive Hamiltonian of the particles in first quantized form. 
A quantum mean field approximation is then formulated 
in section III in order to estimate the transition point and 
asses the importance of quantum fiuctuations within the 
interaction range relevant to ultra-cold dipolar molecules 
or Rydberg atoms. In section IV we derive an effective 
long-wavelength theory of the Zigzag transition starting 
from the microscopic Hamiltonian of the dipolar chain. 
Using the long- wave description we discuss the nature of 
the critical point and the non local order parameter asso- 
ciated with the zigzag phase. The low energy edge states 
imposed by the non-local order in the bulk are discussed 
in section V. In particular we predict direct signatures 
of the zero energy edge states in local tunneling spectra. 
In section VI we demonstrate the key properties of the 
zigzag chain, in particular the presence of string order, in 
a numerical calculation using the density matrix renor- 
malization group (DMRG) method. Finally section VII 
provides a summary and discussion of the results. 



II. MODEL HAMILTONIAN 



We consider a dipolar quantum gas tightly confined in 
a tube-shaped trapping potential. The particles of the 
gas may be either fermionic or Bosonic and we assume 
that their dipoles are polarized by a large electric field (or 
magnetic field in the case of atoms with a large magnetic 
dipole moment) in the direction z perpendicular to the 



tube axis. Such a setup is described by the Hamiltonian 
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where P^, and Vi are the first quantized momenta, 
position and transverse radial coordinates of the particles 
in the trap, uj^ is the transverse harmonic confinement 
frequency m is the mass and d the dipole moment of the 
particles. 

A simple way to analyze the Hamiltonian (1) is through 
a classical mean field theory. This approach consists in 
minimizing the potential energy while neglecting the ki- 
netic terms that induce fiuctuation in the positions. A 
zigzag transition is captured by a Landau-like expan- 
sion of the configuration energy in powers of the stag- 
gered distortion r of the particles.^ For transverse con- 
fi nement frequency uj^ below the critical frequency uJc — 

279 C (5) d'^p^f 8m, the classical energy is minimized by a 
non vanishing distortion r = Hi pQ^^/l — {oo^/oJcY, where 
n = V'186C(5)/3175C(7) ^ 0.245, C(^) is the Riemann 
zeta function, and po is the particle density. In real- 
ity, this criterion only marks the characteristic frequency 
around which a local zigzag distortion begins to develop. 
Quantum fiuctuations allow the zigzag to twist and turn 
and may ultimately destroy the long range order. 

To describe the soft fiuctuations of the zigzag we as- 
sume that each particle is distorted off the axis of the 
trap by the fixed classical value r, but is free to rotate 
around the axis by the angle ip. All our expansions below 
are valid when the distortion r is small compared to the 
inter-particle distance, that is f = rpo should be taken 
to be a small dimensionless parameter. In section (HI) 
we will see that in the relevant range of parameters for 
dipolar molecules this is always satisfied near the zigzag 
quantum phase transition. 

The effective Hamiltonian which describes the linear 
motion of particles along the trap axis, rotations around 
the axis at a fixed radius r, and dipolar interactions be- 
tween particles is given by 



J>1 
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Jij i cos{(pi -(pj)^^ cos{(pi ^(fj) cos2(/?, 
j>i ^ 

Here Pi and Xi are the linear momenta and position oper- 
ators and Li and (pi are the angular momenta (along the 
axis) and its conjugate angle operator of the z'th par- 
ticle. The first two terms in the Hamiltonian (2) de- 
scribe the particle motion and interactions along the trap 
axis, which give rise to density fiuctuations. The rest of 
the Hamiltonian describes the fiuctuations that drive the 
Ising transition. Apart from the last term, the angular 
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part of the Hamiltonian looks like a quantum xy model. 
The dipolar interaction, through the cosine coupling fa- 
vors a staggered arrangement of the molecules off the 
axis, while the angular kinetic energy delocalizes the an- 
gle and thereby acts to disorder the zigzag. The last term 
breaks the U(l) angle symmetry down to Z2 by favoring 
the up and down (0 and tt) distortions of the particles. 
If the trap potential has cylindrical symmetry then the 
preferred axis is set only by the external electric field 
which polarizes the dipoles. In this case we have v = 1. 
It is possible to change v by tuning the ratio between the 
transverse trap frequencies parallel and perpendicular to 
the electric field. 

The other coupling constants appearing in the effec- 
tive Hamiltonian (2) are calculated directly from the full 
Hamiltonian (1). They are given by Jij = J/{po\xi—Xj\)^ 
where J = GcPr'^pQ and / = mr^. Note that the strength 
of the coupling Jij depends on the particles positions 
such that it produces a coupling between the angular and 
longitudinal degrees of freedom. At this point it is con- 
venient to express J and / in terms of two important di- 
mensionless parameters that can be independently tuned 
in the system: (i) r = rpo is the ratio of the classical 
distortion radius to the average inter particle distance 
and (ii) The dimensionless dipolar interaction strength 
Rs = d?pQm/h^ is the ratio of the typical dipolar energy 
Cd = pQ to the typical kinetic energy cq = Ji? Pq/tti. Rg 
can be varied for example by changing the dipole mo- 
ment d using an external electric field, while f is tuned 
by varying the transverse trap frequency uj±. We can 
express the Hamiltonian parameters using these dimen- 
sionless numbers as: J = Gi^gf^eo and I = h'^r'^ /cq. 



III. ESTIMATION OF THE TRANSITION 
POINT 



Within a classical analysis^ the zigzag transition takes 
place at a critical value of the transverse trap frequency 
given by ujc = 6.014:^/R^eo/h. Quantum fluctuations 
driven by the kinetic energy in the Hamiltonian (1) or 
(2) will lead to a transition at a smaller value of the 
transverse trap frequency. With increasing strength of 
interaction Rg the relative importance of the kinetic en- 
ergy decreases and we expect the critical uj± to approach 
the classical value. 

We shall now formulate a quantum mean field theory 
of the effective Hamiltonian (2) in order to estimate the 
value of the parameters at the transition point and to 
asses the importance of quantum fluctuations. For this 
purpose we freeze the axial degrees of freedom and focus 
on the angular part of the Hamiltonian (2). In addition 
since Jij decays rapidly with distance we consider only 
nearest neighbor interactions. 

The mean field approximation consists of decoupling 
the interaction term in the Hamiltonian (2) to get the 



local Hamiltonian 
Hmf = ^ 
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This is supplemented by the self consistency condition 



a {cos (pi)MF ' 
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Note that the Hamiltonian (3) is fully quantum mechani- 
cal in the sense that the non-commuting variables (fi and 
Li are both present in it, in the interaction and kinetic 
terms respectively. Furthermore the transition is a result 
of competition between these two terms, and not driven 
by interaction alone as in the classical zigzag instability. 
On the other hand this local mean field scheme neglects 
the spatial structure of the fluctuations and hence can- 
not capture the universal critical behavior close to the 
transition. Below we estimate the fluctuation region in 
which the mean field theory fails. 

Within the mean field approximation, the system de- 
velops a non zero order parameter a at the critical value 
of the couplings IJ = a{i'). For example we find 
a(l) = 0.31, a(0.5) = 0.33 and a(2) = 0.28. While 
larger values of the eccentricity can be achieved in exper- 
iment, our quantum rotor representation is not suited to 
describe them quantitatively. 

Using the relations IJ = dRgV^ and r = 
tvy^l — {oj^/uJcY we infer a phase boundary 
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The phase diagram is presented in Fig. 2 for the case of a 
symmetric trap (z^ = 1). We see that the transition point 
indeed approaches the classical value of uJc in the strong 
interaction limit. Reasonable values for the dimension- 
less interaction Rg expected in real systems can be es- 
timated from the typical dipolar moments ranging from 
d 0.05 D for magnetic moments through d :^ 0.5 D for 
polarized molecules and up to (i :^ 10 D for atoms weakly 
dressed with Rydberg excitations. Assuming densities of 
po ^ lO^cm"^, we get Rg between 0.05 and 130. From 
the phase diagram (Fig. 2) it is clear that in this regime 
quantum fluctuations shift the transition to far be- 
low the classical value Uc. This implies a huge effect of 
quantum fluctuations on the transition in the relevant 
parameter regime. 

We asses the accuracy of the quantum mean field ap- 
proximation used to compute the transition point by ap- 
plying the appropriate Ginzburg criterion (5cr^)/(cr)^ ^ 
1. The shift of the mean field transition line from the 
line on which the Ginzburg criterion is satisfied (black 
short-dashed line in Fig. 2) estimates the error in deter- 
mination of the critical point. In terms of the dimension- 
less coupling a = IJ the Ginzburg criterion is satisfied at 
a ~ 0.4 (for u = 1) compared to the mean field transition 
found at a = 0.31. 
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FIG. 2: Phase diagram of polar molecules in an elon- 
gated trap in the space of dimensionless interaction constant, 
Rs = dPmpQ and transverse confinement uj±/uJc^ is the 
critical transverse frequency in the classical limit. The solid 
curve marks the zigzag transition (5) computed within the 
local mean field approximation for v — 1 i.e. for isotropic 
transverse confinement. The dashed line corresponds to the 
Ginzburg criterion below which the we expect the mean field 
approximation to break down. Realistic values of Rs range 
between 0.05 and 130 as marked by the grey band. Note that 
the quantum disordered zigzag and the linear chain are not 
sharply distinct. There is a smooth crossover between the two 
regimes. 



IV. LONG WAVELENGTH THEORY 

To capture the universal properties of the zigzag tran- 
sition and the effects of the low energy longitudinal fluc- 
tuations, we derive an effective long wavelength descrip- 
tion of the Hamiltonian (2). The longitudinal modes are 
treated in a standard way. We replace the particle po- 
sition Xj with the smooth displacement field (i)p{xj) and 
the particle momenta by the conjugate field ^dxOp{xj). 
These replacements allow to take the continuum limit of 
the first line in (2) to obtain 



hup 



dx 



Kp{djpf 



K. 
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where Kp = 7t/\^12Rs and Up = \^12R^{hpo/m). From 
here on we rescale the coordinate x by the inter-particle 
distance p^^ and rescale energies by the typical kinetic 
energy cq = h^pQ/m. After this transformation the ve- 
locity becomes dimensionless Up = ^/12Rs. 

The long wavelength limit of the angular Hamiltonian 
should be taken with more care. Because the angle has 
a staggered arrangement it cannot be directly replaced 
by a continuum field (f>{xi). Instead we must make a 
staggered transformation 



{ipi + T^Ni) mod 27r , 



(7) 



where A^^ = Q{xi — Xj) counts the number of particles 
to the left of the z'th particle. Thus the new variable is 



slowly varying in space, with the price of being a non lo- 
cal operator. Now we take the continuum limit through 
the substitution (j)^{xj) and Lj :^dxO(j{xj). 

In addition, noting that Ni irpox — (j)p{x) we can ex- 
press the original local Ising distortion as a composite 
of the two slowly varying fields = r cos{(p{x)) = 

rcos{(j)^{x) + (j)p{x) - irpox). 

The long wavelength Hamiltonian in the angular sec- 
tor, including only the most relevant terms is given by 



Ho 



T~Lc 



27r 



-{dxC^af 



-A 



jdx K,{dJ,f 
j dx{—gi cos 2(j)cr + g2 cos 20cr) , 
dx dx(t)p (1 — cos 2(j)a) ' 



(8) 



(9) 



Here again we have rescaled the coordinate x ^ xpo and 
energies e e/eo in order to work with dimensionless 
coefficients. The term gi originates directly from the last 
term in (2), from which its (bare) value can be estimated 
to be ^1 = J = GRgv'^. The term g2 accounts for the com- 
pactness of the original angle variable. The value of the 
coefficient is the fugacity of kink-anti-kink pairs in (per- 
It is estimated from the action associated with tunneling 
between the two potential minima 0cr = and tt through 
the barrier Jcos2^cr-^^ Within the WKB approximation 



we obtain g2 



/ij 



'exp{-8VRsr^)' In ad- 



dition we extract from (2) the bare value of the Luttinger 
parameter = tta 

ip/V2 = V6R~s 



^TJ = 7r\^6r^Rs and the velocity 

Finally the term (9) describes the coupling between 
the longitudinal phonons and the zigzag degrees of free- 
dom. It is generated because the energy gain associated 
with a zigzag arrangement, proportional to J^j, is mod- 
ulated with the distance Xi — Xj between neighboring 
atoms along the chain. The field dxc/) which couples to 
the zigzag energy in (9) is simply the long wavelength 



limit of 



The value of the coupling constant is 



again extracted directly from (2), A = Ibir Rsv'^. Other 
coupling terms between the two sectors have been omit- 
ted since they turn out to have lower scaling dimension 
at the critical point which describes the zigzag transition. 

We note that our derivation of the field theory was car- 
ried out in the strongly interacting limit ^ 1, where 
the particles almost form a Wigner crystal. In this case 
the low energy Hamiltonian is essentially independent of 
the statistics (Bosonic or fermionic) of the constituent 
particles. As we shall see later, the differences due to 
particle statistics arise when we try to probe the system 
by inserting or extracting single particles to it. 

Having derived the effective field theory we are now 
in position to discuss the resulting quantum phases and 
phase transitions. Let us start with the phases. The 
Hamiltonian Ha of the angular degrees of freedom de- 
scribes a competition between two cosine terms in a 
regime where they are both relevant with respect to the 
quadratic part of the Hamiltonian. In the zigzag ordered 
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phase the term gi is dominant. The field is then 
pinned to one of the two minima at or tt, and the Ising- 
like order parameter {cos (per (x)) takes a non vanishing 
value +1 or —1. On the other hand when g2 becomes 
dominant the dual variable 6a- is pinned, leading to pro- 
liferation of phase slips in (j)s and disordering of the Ising 
field. 

The representation of the order parameter in terms of 
the Bosonic field (j)^ conceals the fact that the order in 
the zigzag phase is non-local. It consists of an infinite 
string operator when written in terms of the actual par- 
ticle distortion = rcos(p{x) 



(rcos0^(x)) = (e 



_ /J^ I-^ p(x) dx' 



(10) 



where p{x) is the particle density. The string order pa- 
rameter captures the fact that the phase of the zigzag, 
namely an up or down distortion, is alternating in the 
particle frame but is not fixed to a position in the lab 
frame. The direct correlation function of the transverse 
distortions cf{x) = r cos (f>{x) decays as a power-law along 
the chain due to disordering of the positional order: 



(a(x)(j(0)) ^ cos(7rpo^) (cosher) \x\ 



(11) 



Let us move on to discuss the critical point separat- 
ing the zigzag from the disordered phase. If we ignore 
for the moment the coupling A between the zigzag de- 
grees of freedom and the longitudinal phonons, then by 
symmetry considerations the effective Hamiltonian (8) 
is expected to describe an Ising critical point. When 
Kfj = 1 the problem becomes particularly simple, since 
we can immediately express l-L^ as a model of two in- 
dependent Majorana Fermions or alternatively a single 
Dirac Fermion:^^'^^ 



e". (12) 



Here Ua- = u^^ = ^l)^ ^ Majorana spinor in 
the chiral basis, and r^'^ are Pauli matrices that act in 
this basis. The masses of the two Majorana modes are 
given by A± = 7r(^2 ^ The transition at gi = g2 is 
then described by a single massless Majorana mode 
while the other gapped mode can be ignored. This is 
precisely the critical theory of the transverse field Ising 
model. The relation between the low energy Majorana 
mode and the Bosonic fields in the low energy limit is 



\/po/ir ■■ cos {9a + (per) ■ 



(13) 
(14) 



From the mean field analysis presented in section (III) 
we expect the transition to occur at a bare value of 
close to, but not necessarily precisely one. A deviation 
from K(j — 1 (see appendix A) gives rise to an interaction 
that couples the two sectors of Majorana Fermions 



(15) 



where V = 27rUa- — K^). In addition there is an 

increase in the velocity of the Majorana modes, which 
now becomes Ua = {^a^ + ^a)- 

Because the interaction couples the critical modes to 
a gapped excitation it is perturbatively irrelevant at the 
critical point. In other words the double sine- Gordon 
model (8) flows under renormalization to the self dual 
(Ising) point = 1. The residual effect of the interac- 
tion is to shift the transition away from the point gi = g2- 
This shift can be estimated using a mean field decoupling 
ofC^CZ fromC^?+ in (15). 

So far we have ignored the coupling A between the Ising 
degrees of freedom and the longitudinal phonons. The 
effect of this coupling, which a priori breaks the confor- 
mal invariance and Lorentz-symmetry of the Ising criti- 
cal point, was addressed in Ref. [12] using a perturbative 
renormalization group analysis in A. To one loop order, 
A was found to be irrelevant if the sound velocity Up is 
initially larger than the velocity of the Majorana mode 
Ucr. The slow flow toward A = has a peculiar effect 
on the fixed point. The ratio between the two veloci- 
ties U(j/up approaches unity, thereby restoring Lorentz 
symmetry at the fixed point. However, at the same time 
the velocities themselves flow to zero. In the opposite 
regime, Ua > Up the coupling A is relevant. The nature 
of the transition in this case is not well understood and 
furthermore corrections beyond one loop can change the 
picture. We would like to point out however that a per- 
turbative expansion in A is justified for the bare values 
of the parameters dictated by our system. The small pa- 
rameter of the expansion is \^Kp/ {u^Up), which in our 
case scales as l/^/R^ <C 1. 

An experiment can in principle explore both of the 
above regimes by tuning the ratio Ua/up. We have shown 
in section HI that changing the eccentricity of the trans- 
verse confinement leads to change of the critical value of 
RgV^ and therefore also of the bare Luttinger parame- 
ter Kfj near the transition point. For = 1 the ratio 
Ua/up = 1/a/2 < 1, while by changing (bare) the 
ratio can be tuned above 1. 

Finally, it is interesting to note the difference between 
our model of the zigzag transition and the problem of a 
two leg Fermion ladder discussed in Ref. [14]. As we have 
mentioned above, such a Fermion model should not be 
essentially different from ours in the strong interaction 
limit. However Ref. [14] assumes the weak coupling per- 
spective, taking the Fermi energy to be much larger than 
both the interaction and the splitting between the two 
sub-bands. This leads after Bosonization (using our con- 
ventions) to a double sine-Gordon model similar to (8) 
in the "spin" (antisymmetric) sector, but with the term 
g2 cos 26 replaced by a cos 46 term, which corresponds to 
Cooper pairing. Accordingly the phase competing with 
the zigzag in the weak coupling limit is a Superconduc- 
tor (i.e. quasi long range order in the pairing field), while 
in our case the competing phase is unpaired. The criti- 
cal theory which describes the transition into the paired 
phase is different from the transition to the the unpaired 
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FIG. 3: States of dipolar particles along an elongated trap with soft harmonic confinement. The inhomogeneous density along 
the trap acts as a tuning parameter of the zigzag transition. Here the particles in the dense middle section are in the zigzag 
phase while the wings of the cloud are in the disordered phase. Majorana zero modes are exponentially localized near the 
interfaces between the two phases. 



state found in our case. Note in particular that the self 
dual point of the double sine-Gordon model with the pair- 
ing term is at K^r = 2 rather than K(j = 1 in our case. 



V. ZERO ENERGY EDGE STATES 

The non-local order has interesting consequences for 
the physics of the edge. Let us consider a physical model 
of the edge as an interface between the string-ordered 
zigzag phase and the disordered phase, created by a spa- 
tially dependent gap parameter A{x) that changes sign 
on the interface. This is in fact a realistic model for dipo- 
lar particles in a harmonic trap (see FIG. 3). The den- 
sity which slowly decreases away from the center along 
the trap is a tuning parameter for the zigzag instabil- 
ity. Therefore, the dense middle section may be in the 
zigzag phase, the wings of the cloud are in the disordered 
phase, and there is necessarily an interface between the 
two regions. 

The Hamiltonian (12) is formally identical to a one 
dimensional Superconductor of spinless Fermions. This 
model has a zero-energy Majorana mode localized at 
the interface, where the gap function changes sign.^^ 
The zero-mode has a particularly simple structure if A 
changes abruptly across the interface from Aq to — Aq 
(see appendix B). It is then given by the Majorana oper- 
ator 



dxe 2z ^ (r{x)) 



(16) 



where / = 'Ucr/4Ao is the correlation length. Bulk states 
of (12) are above the gap Aq. 

Here we should remark on the difference between the 
Majorana edge modes of the zigzag phase and those found 
in two other one-dimensional systems. Consider first the 
transverse field Ising model. The representation of this 
model in terms of Jordan- Wigner Fermions shows Majo- 
rana zero modes localized at the two edges. But when 
mapped back to the original Ising degrees of freedom, 
these two modes simply span the Hilbert space of the 
two broken symmetry states in the bulk. In our case, be- 
cause there is no broken symmetry in the bulk the zero 



modes are truly localized at the edge, as we will show 
explicitly below. 

There is also an important difference between the Ma- 
jorana edge modes discussed here and those of a spinless 
Fermion model in which the Fermions are the local de- 
grees of freedom. In our case, the Majorana modes ap- 
pear only in the spin-sector, where they describe a double 
degeneracy of the spectrum for a given charge parity of 
the zigzag. For even parity, the two states can be under- 
stood as either a zigzag that starts with an up distortion 
and ends down, or the opposite configuration. For odd 
parity the two states correspond to a zigzag with an up or 
down distortion at both ends. However in a system with 
an interface, where the parity is not well defined, all four 
states should exist and may be combined into completely 
local excitations. This is in contrast to the non-local na- 
ture of the edge modes in the Fermionic wire.^^ 

We shall look for physical signatures of the edge states 
in the local tunneling density of states (DOS).^^ Consider 
a probe that can tunnel particles into or out of the system 
at a desired point x and energy uj. Such local probes have 
been proposed^^, and are now becoming a reality in cold 
atomic systems. ^^'^^ We also assume that the position of 
the tunneling tip is asymmetric with respect to the tube 
axis so that it can insert a particle either displaced up or 
{s = t) down {s = D in the tube. The tunneling rate is 
proportional to the local spectral function 



As{uj,x) 



Im 



(17) 



where the ± subscript denotes anti-commutator and com- 
mutator for Fermions and Bosons respectively. To find 
the low energy behavior of the spectral function we seek 
the representation of the particle field operators in terms 
of the Bosonic density modes (j)p and Op and the low en- 
ergy Majorana Fermion ^. This will depend in a crucial 
way on the species, i.e. Bosonic or Fermionic, of the dipo- 
lar particles. In what follows we separately treat each of 
these cases. 

Bosonic particles- Let us start from the standard 
Bosonization identity for Bosonic particles in two chains 
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labeled by s = t (-) 

OO 



m= — oo 



We now canonically transform to a representation 
in terms of symmetric ("charge") and anti-symmetric 
("spin") fields through Os = Op -\- sOa- and (j)s = {(l)p + 
S(\)(j)l2. Keeping the most relevant terms m = 0, ±1 
and using the Bosonized form (13,14) of the Majorana 
operators we obtain 



(19) 



where /3o and /3i are non- universal coefficients. The 
first row in the operator (19) corresponds to the long- 
wavelength components of the Boson. The second row 
describes the components with wavelength comparable 
to the inter-particle spacing. In both terms insertion of a 
Boson must cause a disruption in the zigzag order. In the 
first case the operator e^^^'^ creates a kink in the Ising 
order-parameter, while in the second case the disruption 
is affected by the Majorana operator. 

Since we are interested in the low energy contributions 
near the interface we can further simplify the Bose oper- 
ator. First the factor e^^ can be replaced by its expec- 
tation value, which decays exponentially with distance 
from the interface in the ordered side (recall that (e^ ) is 
the "dual Ising order-paramete") . Second, the Majorana 
operators ^.l.r can be replaced by the contribution to 

_ I I 

them from the zero energy edge mode, ~ e ^ '^'70- 
From here we can easily compute the local spectral func- 
tion at low energies (see appendix. D) 



A{E,x) 



1 



i x/l 



Bie-I^l/'IE^I^" ) \E 



(20) 

where 5o,i are non-universal constants. The leading 
(first) term originates from the long wavelength contri- 
bution to the Bose operator. Only the sub-leading term 
Bi stems from the Majorana zero mode. 

Fermionic particles - Unlike with Bosonic particles, 
we cannot access the zigzag phase of Fermions starting 
from weakly coupled Luttinger liquids. A better start- 
ing point to obtain the low-energy limit of the Fermion 
field is a two band model, including the ground state and 
first transverse excitation in the tube.^^^^ Details of the 
formulation are left to appendix C, however the end re- 
sult is simple. A Fermion inserted into the lower band is 
described by 



4 



"0 VPo 



E 



(21) 



This is directly analogous to the first term in (19), except 
that the Fermion is not inserted at zero momentum but 
rather at the Fermi points of the Luttinger liquid that 



describes the breathing modes of the zigzag. Insertion 
of a Fermion must be accompanied by creation of a kink 
in the zigzag order, which is implemented by the factor 
g±i Because on both sides of the transition there is a 
gap to a second band, the Fermion can also be inserted 
at zero momentum above the gap. This "second band" 
Fermion is directly proportional to the Majorana modes 
(see appendix C) 



(22) 



A Fermion distorted up or down in the tube is now de- 
scribed by a superposition of the symmetric and antisym- 
metric transverse states so that '^f,s = (V^o + 5?/^i)/>/2. 

With the Fermion operators at hand we are now in 
position to compute the local spectral function related 
to the tunneling rate of up or down distorted Fermions 
into a point x. The result is again a power-law 



A{E,x) 



Ao 



1 



■ Ai e 



-\x\/l 



\E\ 



(23) 

As in the Bosonic case, the first term describes the "pene- 
tration" of the disorder-parameter into the ordered phase 
near the interface. The second term is the contribution of 
the Majorana mode. Note however that contrary to the 
Bosonic case, here the Majorana mode gives the domi- 
nant contribution to the low energy tunneling rate near 
the interface. It is interesting to note that Eq. (23) shows 
the same leading behavior with E as found in Ref. [14] 
for tunneling into the edge of the zigzag phase in the 
weak coupling limit. This is in spite of the very different 
modeling and physical picture of the edge in the two sce- 
narios. From our analysis it is clear that this behavior 
stems from a robust zero energy state with topological 
origin. 



VI. NUMERICAL RESULTS 

We now turn to a numerical calculation of a concrete 
microscopic model using DMRG^^ and verify that the 
long range correlations are indeed well described by the 
field theoretical analysis presented above. For this pur- 
pose we consider a lattice model which describes similar 
physics as the original Hamiltonian (2). Specifically we 
take a model of hard-core Bosons on a two leg ladder at 
low incommensurate filling (see Fig. 4. (a)): 



HLadder = ^ - t|| ^ (4,a^i+l,a + h.C.) (24) 
i=l ^ a=l,2 



a=l,2 



where b] ^ create a (hard-core) Boson on site i of leg a = 
1,2. The two legs of the ladder represent the preferred 
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rja r/a 



FIG. 4: (Color online) DMRG calculation of correlations in 
the two- leg ladder model (24). The two legs of the ladder 
illustrated in panel (a), correspond to the up and down dis- 
tortion. V± is taken to be the largest energy scale to exclude 
double occupancy of a rung, and the zigzag transition is tuned 
by varying the ratio t±/V\\. (b) The component with spatial 
frequency tt/Oo of the Ising (red dashed) and string (blue solid) 
correlation functions in the ordered phase on a Log-Log plot. 
The string correlations display true long range order while the 
Ising correlations decay as a power-law. (c) The same correla- 
tion functions calculated in the disordered phase plotted on a 
Log-Linear plot, where both are seen to decay exponentially. 

up and down distortion of a particle. The corresponding 
Ising variable is the relative density cfi = n^^i —ni^2- The 
repulsive interaction Vj| between particles on the same 
leg favors a zigzag arrangement to minimize the interac- 
tion energy, just as the dipolar interaction does in the 
original system. At the same time the hopping on a rung 
t± counters that ordering tendency. The latter is di- 
rectly related to the tunneling matrix element between 
the two favored states of the rotor = and tt. Tuning 
t± to drive the phase transition is analogous to varying 
the transverse trap frequency in the original model. The 
hopping t|| along the ladder drives the longitudinal quan- 
tum fluctuations. At incommensurate filling, just as in 
the continuum, these fluctuations prevent crystalline or- 
der from forming. Finally we fix a large repulsion V± 
(^ V\\^t\\^tj_) to suppress double occupation of a rung. 

To characterize the ground state we compute the stag- 
gered Ising correlation function and string correlation 
function, which are defined respectively as: 

Clin - Vj) = cos {7rpo{ri - rj)) {ai&j) 

Cstr{n-Tj) = {a,e'^^^=o^^a^). (25) 

Results of two calculations which differ only by the value 
of the tuning parameter are presented in Fig. 4. The 
lattice size in both cases is L = 256 unit cells (rungs), 
the filling is po ^ 1/3 per unit cell, and V\\ = t\\ = 1. 
The transverse hopping is chosen to be t± = 0.05 in 



one calculation (panel (b) of the figure) and t± = 0.8 
in the other (panel (c)), so that the system is in the 
zigzag ordered phase in the first and in the disordered 
phase in the second calculation. As anticipated, the Ising 
correlations decay like a power law of distance in the 
ordered phase while the string correlation saturates to a 
finite value. When the transverse tunneling is increased 
to t± = 0.8 both correlations decay exponentially with 
distance. 



VII. CONCLUSION 

The zigzag transition observed in chains of trapped 
ions^ is a mechanical distortion that occurs when the ion 
chain is compressed or the transverse confinement de- 
creased beyond a critical point. The transition is well 
described by a classical theory neglecting all quantum 
(as well as thermal) fluctuations.^ In this paper we have 
shown how the analogous transition takes place in an 
elongated trap of ultracold dipolar particles, where quan- 
tum fluctuations are greatly enhanced compared to the 
ion system. 

Quantum fluctuations are found to have a profound ef- 
fect on the transition and on the nature of the ordered 
phase. In particular, the Z2 symmetry of the system, 
which is broken in the classical zigzag phase, is restored 
by the quantum fluctuations of particle positions along 
the trap. Instead, the quantum zigzag phase is char- 
acterized by a non-local, Ising string order parameter. 
We also find zero energy edge states, localized at inter- 
faces between the zigzag phase and the disordered (lin- 
ear) phase, that are a direct consequence of the non local 
order in the bulk. It is notable that such non trivial cor- 
relations are found in the most natural setting for polar 
molecules with no need for special engineering of inter- 
actions. 

We have calculated the experimental signatures of the 
edge states that are expected to be seen in tunneling ex- 
periments. It would also be interesting to observe the 
bulk string order directly using in situ detection of the 
particle positions. From the fundamental perspective, 
our analysis provides another example of a rather rare 
class of systems which exhibit topological order (or at 
least non-local string order) in a gapless phase. Other 
examples in this class include the Haldane phase of spin- 
3/2 chains and a Haldane liquid phase predicted to oc- 
cur for two component dipolar Fermions in an optical 
lattice.2^ The zigzag phase is somewhat special in this 
class in that it occurs in presence of full (continuous) 
translational symmetry rather than in a lattice. 
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Appendix B: The Majorana edge state 



Appendix A: Re-Fermionization of the double 
sine-Gordon model 

Here we review the re-Fermionization scheme for the 
double sine-Gordon model^^'^^'^^(8). Let us start by re- 
writing the Hamiltonian (8) as follows: 

= 1^ J dx [{d,e,f + {d,<p,f] (Ai) 



2n 2K, 



dx{gi cos 20cr + 92 cos 20^) • 



Using the identities 



^/^.^-HeM .^^-+^+^ (A3) 

we find that the re-Fermionized version of the first and 
third terms in this Hamiltonian (Al) give the quadratic 
Majorana model 



Hi") 



, (A4) 



where ^ 



± _ 



= + Kl)/2K, and 
A± = 7r((/2 ±5i)/Po- The re-Fermionized version of the 
second term in (Al) has the form of an interaction 

H^^'^ = V j dxCR^L^tCL- (A5) 

where V = 2t:U(j{1 — KD/2K^. As explained in section 
IV this interaction is irrelevant in the renormalization 
group sense and therefore not expected to alter the criti- 
cal properties. It can however affect a shift of the critical 
point from the naive value A_ =0. As explained in the 
text this shift can be estimated by decoupling the inter- 
action term V so that the renormalized values of and 
A_ become mean field parameters to be determined self 
consistently. If V is small compared to the shift in 
A_ can be easily obtained perturbatively by substituting 
^R^L interaction term by their expectation value 

in the unperturbed ground state 



log 




Then the renormalized value ofA_ isA_ = A_ — 
i(^^^^)V. In the limit A+ ^ Po^cr this approaches 
A_ = A_ + V. 



In this appendix we solve the Bogoliubov-deGennes 
(BdG) equations in the presence of an edge and show 
that there is a single zero energy solution. For simplicity 
we take a sharp change in the mass term 

..xJAq x<0=^ disorderd . 
^-^^^~\-Ao x>0 ^ orderd * ^^^^ 

When Kfj = 1 the critical point of the Ising model (8) is 
at A_ = 0, therefore, equation (Bl) defines a boundary 
between an ordered phase for x > and a disordered one 
for X < 0. The BdG equations for this situation are given 

by 



-i Y^.T- + A_(x)t^J xe{x) = Exe{x) (B2) 



where Xe{x) = {ue{x)^ve{x)}^ is the eigenstate of the 
equation (B2) with energy E . The physical solutions of 
(B2) are at energies above the gap (£' > Aq) except for 
the single zero energy solution 



Xo(^) 



2Ao 



-1 



(B3) 



The solution of the BdG equations immediately gives us 
the quasi-particle operators 



= J dx [ue{x)^r{x) - ve{x)£,l{x)] 
7e = J dx [ue{x) ^r{x) - v%{x) ^l{x)] 
and especially the localized Majorana operator 



70 = 76 



2Ao 



(B4) 
(B5) 

(B6) 



Appendix C: Low energy limit for Fermionic dipoles 

In case of Fermions it is simpler to derive the Bosoniza- 
tion identities starting from the weak coupling picture 
considered in Refs. [10-12]. Here we will review the 
derivation of the low energy theory in this case and draw 
analogies to the field theory derived above (8) starting 
from strongly coupled dipolar particles. Guided by this 
low-energy theory we will give a phenomenological ex- 
pression for the field operators. 

In Ref. [12] the two lowest sub-bands of the transverse 
confining potential are considered. The zigzag transition 
takes place when the chemical potential is placed at the 
bottom of the second sub-band. At this point the group 
velocity of the second sub-band is zero and the Bosoniza- 
tion procedure can not be carried out. However, we can 
still Bosonize the first sub-band. The Hamiltonian of the 
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two sub-bands can be written as follows 



H 



27r 



dx 



K{dJo) 



dx ■i/'l 



'2m 



1 

K 



(CI) 



/ 



dx 



-d^^oi^liPi + -^ie'''°iPAi^i + h.c.) 



where S is the difference between the energy of the second 
sub-band and the chemical potential, S = ei — jj.. The 
field operator ipi belongs to the second sub-band and the 
first sub-band has been Bosonized 



• •. • 



(a) 



m • IB • 



TT 



X 



EUr lim 
nr\ — Vor 



P0_ i rkFX -\ (r0o-0o) 



(C2) 



r=R,L 



where is the appropriate Klein factor. The density- 
density interaction g and the pairing u in the Hamiltonian 
(CI) originate in the coulomb interaction.^^ 

In order to decouple the pairing term in the Hamil- 
tonian (CI) the authors of Ref. [12] apply the unitary 
transformation 

jj — f dx' Ooiplipi ^ (C3) 

such that the transformed operators are given by 

^0 = U'^^l^oU - (C4) 



x^-i (r(/)o- 



o)p-^^ f dx'sign{x-x')iplipi 



r=R,L 



and 



The Hamiltonian (CI) then assumes the form 



(C5) 



(C6) 



• ^1 z 



25Ty 



Here the curvature of the second sub-band has been 
neglected compared to the linear pairing term u. We 
have also switched the notations here such that dx(t>p = 
9x4*^ + is the fluctuation in the total density, and 

9p = Oq is the conjugate phase. We would also like to 
note that one can recover the Majorana theory (12) by 
decomposing into it's real and imaginary parts as fol- 
lows 



FIG. 5: (a) and (b) A chain with a soliton that is created by 
a defect in the s =t and s =| state respectively, (c) In both 
cases the field ^o- is shifted by — tt due to the soliton. 



single non-chiral massive Majorana field. Therefore, we 
identify our theory as the strong coupling limit of this 



Apparently, the expressions for the single particle oper- 
ators (C4) and (C5) are sufficient to express the Fermion 
operators in the s =t, i state in the zigzag. However, the 
operator (C4) is written in an inconveniently non-local 
form. We can overcome this inconvenience by replacing 
(C4) with a phenomenological expression. As shown in 
Fig. 5 the insertion of a particle both into the s =t and 
s =1 states in the zigzag shifts the field (j)a- by ±7r. The 
operator that shifts (j)a- by ±7r is e^^^'^. Therefore, we 
replace the non local string with Oa- 



E 

-=R,L 



(C8) 



This operator creates a soliton in the a sector in addition 
to a plasmonic excitation in the charge sector. Similarly, 
the operator has the form of a spinon multiplied by 
the factor e^^p that inserts a charge at zero momentum. 

We can now write expressions for second quantized 
Fermion operators inserted to the up and down positions 
in the tube s ="[,-1 as the appropriate superpositions of 
Fermions from the symmetric and anti symmetric sub- 
bands 



i^) ■ 



(C7) 



The Hamiltonian (C6) and (6), (8) and (9) are identical. 
They both consist of two sectors a Luttinger liquid and 



■S1pi). 



(C9) 
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Appendix D: Tunneling DOS on the edge 
1. Bosons 

The low energy limit of the second quantized Bose op- 
erator is (19) 



27rAo 



l\x\ 



X 

(Dl) 



IS e 



-i((f)p-\-27Tpox) 



70 • 



The Bose operator is used to compute the single par- 
ticle Green's function in the effective field theory (6,8). 
The contribution of the first term in the Bose operator is 
found by noting that the disorder parameter e^*^^ has an 
expectation value in the disordered side of the interface 
that penetrates to the ordered side and decays there as 
^-x/i _ {^^/4/\q ig the correlation length in the ordered 
phase). We can therefore replace the operator e^*^^ by 
this expectation value. The contribution from the sec- 
ond component of the Bose operator is computed using 
the form of the zero-mode (B3) and the Luttinger liq- 
uid which describes the density modes. Putting all this 
together we have 



GB{r,x) = -(T,^B,.(r,x)^;[,^^(0,x)) 



(D2) 



l/^olVo 



1 + 

\x\/l 



2po / 1 + e-^/^ 
1 



1 



PoUp\r\ 



21 \PoUp\r 
where we have used the fact that 

{IV7o(t)7o(0)> = ^Sign(T) . 



(D3) 



To obtain the retarded greens function we analytically 
continue the function (D2) 



GB(t,x) 



(D4) 



sm 



8K, 



1 



21 



l + 2po/ 1 + e-^/^ VPoV. 



1 



-\x\/l 



1 



pQUpt 



The imaginary part of the Fourier transformed version of 
this function gives the local DOS: 



where, 



1 



i x/l 



■Biel^l/'l, 



• (D5) 



Bo - |/3o|Voi 



tt/2 



(1 + 2po r[l/4i^p] \poUp 



and 



2 1 



21 V[Kp^l/AKp] \poUp 



The leading contribution to the local DOS stems from the 
penetration of the disorder field into the ordered phase. 
The sub- leading contribution originates from the Majo- 
rana edge state. 



2. Fermions 

The low energy limit of the Fermion operator was 
found above and is given by 



F,s 



1 

71 



^0 Yl v^e^'^^'^'^e-^^^^^-^^^e^^^- (D6) 



We use this operator in the effective low energy theory 
(6,8) to obtain the imaginary time Green's function in 
the same way as discussed above for the Bosonic case. 
This gives 



Gf{t,x) ^ -Sign(r)x 



(D7) 



|ao|Vo e 



_vTT2^1 + e-^A \poUp\r\ 



-\x\/l 



PoUp\r\ 



To obtain the low-energy local DOS we analytically con- 
tinue this expression and the imaginary part in fourier 
space 

A^(c.,x)^^^^£^^^+AieN/'|c.|*-^ (D8) 



where. 



Ao ^ 1^0 1 Po 



and 



l + 2po/r[i^^ + l/4i^^] 



I \^P^4Kp 



1 



21 T[l/4Kp] \poUp 



In this case the leading contribution comes from the lo- 
calized Majorana Fermion, and the sub-leading one stems 
from the penetration of the disorder parameter into the 
ordered phase. 
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